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We present a simple procedure to obtain all de Sitter solutions in the ghost-free 
massive gravity theory by using the Gordon ansatz. For these solutions the physi¬ 
cal metric can be conveniently viewed as describing a hyperboloid in 5D Minkowski 
space, while the flat reference metric depends on the Stuckelberg field T[t , r) that sat¬ 
isfies the equation (<9jT) 2 — ( d r T ) 2 = 1. This equation has infinitely many solutions, 
hence there are infinitely many de Sitter vacua with different physical properties. 
Only the simplest solution with T = t has been previously studied since it is mani¬ 
festly homogeneous and isotropic, but it is unstable. However, other solutions could 
be stable. We require the timelike isometry to be common for both metrics, and this 
gives physically distinguished solutions since only for them the canonical energy is 
time-independent. We conjecture that these solutions minimize the energy and are 
therefore stable. We also show that in some cases solutions can be homogeneous and 
isotropic in a non-manifest way such that their symmetries are not obvious. All of 
this suggests that the theory may admit viable cosmologies. 
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I. INTRODUCTION 

The discovery of the ghost-free massive gravity theory by de Rham, Gabadadze, and 
Tolley (dRGT) [1] (see [2] for a review) opens up the possibility to explain the dark energy 
and the cosmic acceleration [3] by a tiny mass of the gravitons. The dRGT field equations 
admit the de Sitter solution with the cosmological constant mimicked by the graviton mass. 
This solution can describe the late time acceleration, but a special analysis is needed to 
decide whether its other properties are physically acceptable. 

A closer look reveals that the de Sitter solution in the dRGT theory is actually not unique, 
and a number of its versions have been found [4] , [5] , [6] , [7] , [8] , [9] . A special attention was 
received by one solution whose physical and reference metrics are of the manifestly homo¬ 
geneous and isotropic Freedman-Lemaitre-Robertson-Walker (FLRW) form [7]. However, a 
detailed analysis revealed that this solution is unstable [10]. For other known solutions only 
the physical metric is manifestly FLRW while the reference metric looks inhomogeneous, 
for which reason they are considered to be less interesting [6]. All of this has reduced the 
interest towards the dRGT theory, the focus shifting towards its extensions, as for example 
the bigravity [11], [12],[13] and other generalizations admitting FLRW solutions [14], 

However, we would like to argue that it may be premature to abandon the dRGT theory 
on the basis of negative evidence obtained from just one solution, because the theory admits 
infinitely many other solutions that could be physically interesting. They all have the same 
physical (de Sitter) metric but different values of the reference metric depending on the 
Stuckelberg field T(t,r ) subject to a complicated differential equation [6], [8], [9]. Below we 
shall describe a simple way to obtain these solutions by applying the Gordon ansatz [15] 
and using the global embedding coordinates. The T-equation then assumes a simple form, 
(d t T) 2 — (d r T) 2 = 1, whose essentially general solution is known. The simplest solution T = t 
[7] is unstable [10] but other solutions could be stable. One can choose T such that both 
metrics are invariant under the timelike isometry, which gives special solutions since only for 
them the canonical energy is time independent. We conjecture that their energy is minimal 
and hence these solutions are stable. We also give explicit examples where the reference 
metric looks inhomogeneous but shares with the physical metric the same translational 
and rotational isometries. Hence, solutions considered to be non-FLRW can actually be 
homogeneous and isotropic. All of this suggests that viable dRGT cosmologies may exist. 
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II. THE DRGT MASSIVE GRAVITY 

The theory is defined on a four-dimensional spacetime manifold endowed with two 
Lorentzian metrics, the physical metric and the flat reference metric = 7]Asd^ A d u ^ B 
with i]ab = diag[—1,1,1,1]. The fields $ A (a;) are sometimes called Stuckelberg scalars. The 
theory is defined by the action 

( 2 . 1 ) 

The metrics and all coordinates are assumed to be dimensionless, the length scale being the 
inverse graviton mass 1/m. The interaction between the two metrics is expressed by a scalar 
function of the tensor 7 /, = a /g^ a f a u where g^ v is the inverse of g /w and the square root is 
understood in the matrix sense, i.e. (y 2 )^ = 7 P Q 7°/ = g^fav If (A = 0,1, 2, 3) are the 
eigenvalues of 7 /, then the interaction potential is given by U = 4 Uk where 4 are 

parameters and Uk are defined by the relations 

Mo = 1, M, = Y, Vi = M, M 2 = Y, AaA b = 1(M 2 - M). 

A A<B 

U 3 = Y. A 7 A B A c = l([ 7 ] 3 -3[ 7 ][ 7 2 ] + 2[ 7 3 ]), 

A<B<C 

M 4 = AoA^Aa = 1([ 7 ] 4 - 6[ 7 ] 2 [7 2 ] + 8[ 7 ][ 7 3 ] + 3[ 7 2 ] 2 - 6[ 7 4 ]). 

Here, using the hat to denote matrices, one defines [ 7 ] = tr(y) = 7 ^ , [y fc ] = tr(y fc ) = ( 7 fe ) M M . 
The parameters 4 can apriori be arbitrary, but if one requires the flat space to be a solution 
of the theory and m to be the Fierz-Pauli mass of the gravitons in flat space, then the five 
bk are expressed in terms of two arbitrary parameters, sometimes called C 3 and C 4 , as 

&o = 4c 3 + c 4 -6, bi = 3 - 3c 3 - c 4 , b 2 = 2c 3 + c 4 - 1, & 3 = -(c 3 + c 4 ), b A = c 4 . (2.2) 

The metric g jlv and the scalars d *" 4 are the variables of the theory. Varying the action (2.1) 
with respect to g, w gives the Einstein equations G IW = T^ v with 

n = {bi + b 2 U 1 + b 3 U 2 + b 4 U:i}Yu ~ + b 3 U 1 + b A U 2 }{ y 2 )4 

+ {b 3 U 0 + & 4 W 4 }( 7 3 )4 - b 4 U 0 ( 7 4 T V -US?. (2.3) 

Varying the action with respect to gives the conservation conditions = 0. These 

are equations for the Stuckelberg scalars, but they are actually not independent and follow 
from the Bianchi identities for the Einstein equations. 
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III. DE SITTER SPACE 

The above field equations admit solutions for which the physical metric is de Sitter. 
Specifically, the de Sitter space can be globally visualized as the hyperboloid 

-a' 0 2 + + y 2 = a 2 ( 3 . 1 ) 

i 

is the 5D Minkowski space with the metric 

ds 2 = -dX 2 + Y dX 2 + dX\. (3.2) 

i 

Rescaling the coordinates, A 0 = at, X t = axi, X 4 = ar with Xi = ( x,y, z), the metric reads 

ds 2 g = a 2 {— dt 2 + dr 2 + dx 2 + dy 2 + dz 2 } 

= a 2 {— dt 2 + dr 2 + dR 2 + R 2 dQ 2 } (3.3) 

where dfl 2 = did 2 + sin 2 Myr and 

R 2 = x 2 + y 2 + z 2 = l + t 2 -r 2 . (3.4) 

Let us choose the flat reference metric as 

ds) = a 2 u 2 {-dT 2 + dX 2 + dY 2 + dZ 2 ) , (3.5) 

where u is a constant and T, A", Y, Z are the Stuckelberg fields. 

It turns out that a, u and T , A, Y, Z can be chosen such that the two metrics fulfill the 

field equations. Specifically, it is sufficient to make sure that they fulfill the following relation 

(the Gordon ansatz) [15], 

= + ( i-C 2 )^K), (3.6) 

where c 0 , ( are functions and 

= -1. (3.7) 

If Eq.(3.6) is fulfilled, then one can see at once that 

7"„ = ZrJZ = »B+(i-c)vx), ( 3 . 8 ) 

since 7 ^ 7 “ = g m fo.v One has ( 7 n Y v = LU n (<5% + (1 — C n )R M f4) and so the energy- 
momentum tensor (2.3) is 

T? = - {PoH - CuP 1 (co)}6^ + u (C - l)P 1 (co)V^V u (3.9) 
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with 

Pm{u) = bm + 2b m+1 u + b m+2 u 2 -, m = 0,1,2. (3.10) 

Let us set u — u where u is a constant chosen such that 


Pi(u) = 0. 


(3.11) 


Then the energy-momentum tensor (3.9) reduces to = — Pq{u)5^ and the Einstein equa¬ 
tions become 

G£ + A<5£ = 0 (3.12) 


with A 


Pq(u). The de Sitter metric (3.3) is a solution of these equations provided that 


J_ _ A _ -Pq(m) 

a 2 ~ 3 “ 3 


(3.13) 


Therefore, the metrics (3.3) and (3.5) will indeed fulfill the field equations if u and a are 
defined by (3.11),(3.13), provided that one can adjust the functions T, X, Y, Z such that the 
Gordon relation (3.6) is fulfilled. 

Let us choose in (3.5) T = T(t,r), X = x, Y = y, Z = z so that the f-metric becomes 
ds 2 = a 2 u 2 {— dT 2 + dx 2 + dy 2 + dz 2 } = a 2 u 2 {— dT 2 + dR 2 + R 2 dVL 2 } . (3-14) 


The two metrics (3.3) and (3.14) are related to each other as 


ds 2 = u 2 ( ds 2 + dt 2 — dr 2 — dT 2 ) . 


This will be compatible with the Gordon relation (3.6) if 


(3.15) 


d„td„t - £>„r£>„r - d„Td„T = (1 - < 2 )VV. 


(3.16) 


Assuming that the indices /i, v correspond to (t, r, $, <p) yields 14 = = 0 and 

(d t T) 2 - 1 = (c 2 - m 2 , 

(d r T) 2 + 1 = (C 2 - 1)V 2 , 

d t Td r T = (C 2 - l)V t V r . (3.17) 


From the first two of these relations one obtains 

(d t T) 2 - 1 


^t 2 = 


C 2 


17 


(■ d r T) 2 + 1 

C 2 -i 


(3.18) 
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while the normalization condition (3.7) determines Finally, inserting (3.18) to the third 
relation in (3.17) yields 


(d t T) 2 (d r T) 2 = ((d t T) 2 - 1 )((d r T) 2 + 1) (3.19) 

and therefore 

(< d t Tf - (d r T) 2 = 1. (3.20) 

This completes the procedure, because V) t and ( are determined by the above formulas and 
the Gordon relation is fulfilled. 

Summarizing, the de Sitter solution in the theory is described by Eqs.(3.3),(3.14) where 
u,a are defined by (3.11),(3.13) and T is a solution of the differential equation (3.20). Since 
there are infinitely many T’s subject to (3.20), there are infinitely many de Sitter solutions. 
They all have the same physical metric (3.3) but differ one from the other by the choice of 
T in the reference metric (3.14). The physical properties of solutions with different T’s, as 
for example their stability, can be different. 

These solutions were actually obtained previously [6],[8],[9], but within a different com¬ 
putation scheme yielding the T-equation in a form that gives little hope to solve it (see 
Eq.(A.3) in the Appendix). Our procedure yields its equivalent form (3.20), which is simple 
and admits a general solution. In addition, by slightly modifying our procedure, we can 
obtain new things. Specifically, it was assumed in the above derivation that both metrics 
have the same spatial 50(3) symmetry, ffowever, let us rather choose 

ds 2 = a 2 u 2 {— dt 2 + dx 2 + dy 2 + dZ 2 } (3-21) 

with Z = Z(r, z), so that the two metrics share the same 50(1,2) symmetry in the t,x,y 
subspace. Repeating the above analysis one obtains 

(< %Z ) 2 + ( d z Zf = 1, (3.22) 

and this gives new solutions. When expressed in the standard spherical coordinates, their 
f-metric will not look spherically symmetric, since for generic Z it has no common with the 
g-metric 50(3) symmetry, although it has its own 50(3) in the x,y,Z space. Below we 
shall mainly be discussing equation (3.20) since the analysis of (3.22) is similar. 
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IV. THE SIMPLEST SOLUTION 

Even though there are infinitely many solutions of Eq.(3.20), almost all known dRGT 
cosmologies reported in the literature correspond to just one simplest solution, 

T — t. (4.1) 

A slightly more general choice is 

T = cosh(£) t + sinh(£) r (4.2) 

with a constant £. However, the value of £ can be changed by a boost in the t,r plane of 
the ambient 5D Minkowski space, which does not affect the g-metric (3.3), hence one can 
set £ = 0 without loss of generality. Rewriting (4.1) in different coordinates gives results 
which look very different, and it has not been recognized that they actually describe the 
same solution. Let us see what happens when this solution is expressed in the standard 
spatially flat, closed, and open coordinate systems. 

1. Flat slicing 

Let us express t,r,R in (3.3) in terms of two new coordinates r and p as 

n 2 P 2 

t = sinh r + ^-e T , r = coshr —— e T , R = e T p . (4.3) 

This solves the constraint (3.4) and transforms the de Sitter metric (3.3) to the standard 
FLRW form with flat spatial sections, 

ds 2 = a 2 {—dr 2 + a 2 (r)(dp 2 + p 2 dfl 2 )}, (4.4) 

where a(r) = e r . The function T = t can be represented as 



This solution was found in Ref. [ 6 ] for given by (2.2) with 03 = 04 = 0, and later for 
arbitrary bk [9] (solution in [ 6 ], [9] contains an integration constant that can be obtained by 
using (4.2) instead of (4.1)). Although the g-metric (4.4) is manifestly homogeneous and 
isotropic, the f-metric (3.14), when expressed in the r, p coordinates, becomes non-diagonal 
and p-dependent, which suggests that it is inhomogeneous. For this reason it is sometimes 
said that the dRGT theory does not admit genuinely homogeneous and isotropic cosmologies 
with flat spatial sections [ 6 ]. However, we shall shortly comment on this. 



2. Closed slicing 


If one chooses in (3.3) 

t = sinh(r), r = cosh(r) cos(p), R = cosh(r) sin(p), (4.6) 

this solves the constraint (3.4) and the de Sitter metric (3.3) assumes the FLRW form with 
closed spatial sections, 

ds 2 g = a 2 {—dr 2 + a 2 (r)(dp 2 + sin 2 (p)dCl 2 )} (4.7) 

with a(r) = cosh(r). These coordinates cover the whole of de Sitter space. The Stuckelberg 
held is T = sinh(r), and the f-metric (3.14) expressed in the r, p coordinates is again non¬ 
diagonal and p-dependent, which suggests that there are no genuinely homogeneous and 
isotropic cosmologies with closed spatial sections either. 

3. Open slicing 


For the open slicing one has 

t — sinh(r) cosh(p), r = cosh(r), R — sinh(r) sinh(p), (4.8) 

and the g-metric becomes 

ds 2 = a 2 {—dr 2 + a 2 {r){dp 2 + skill 2 (p)dQ 2 )} (4.9) 

with a(r) = sinh(r). The Stuckelberg held is T = sinh(r) cosh(p) and the specialty now is 
that the f-metric (3.14) becomes diagonal in the r, p coordinates, 

ds 2 = a 2 u 2 {— cosh (r) 2 dr 2 + a 2 (r){dp 2 + sinh 2 (p)dfl 2 )}. (4.10) 

This solution, discovered in Ref. [7], is broadly viewed as the only homogeneous and isotropic 
dRGT cosmology, because both metrics are manifestly homogeneous and isotropic, so that 
they share the same rotational and translational Killing symmetries. However, this solution 
is completely equivalent to its hat and closed versions. Therefore, the latter also have the 
same common isometries, hence they are all homogeneous and isotropic, although their 
symmetries are not manifest. The conclusion is that sometimes solutions can be FLRW in 
a non-manifest way. 

At the same time, although homogeneous and isotropic, the solution T — t is not static 
whereas the de Sitter space is. Specifically, let us consider the 
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4- Static slicing 


Setting 

t = \Jl — p 2 sinh(r), r = sjl — p 2 cosh(r), R = p (4.11) 

solves the condition (3.4) and reduces the de Sitter metric (3.3) to the static form 

ds 2 g = a 2 Edr 2 + 

with E = 1 — p 2 . The T = t solution then becomes 

T(r,p ) = \/l - p 2 sinh(r), (4.13) 

and it is non-static even in static coordinates. Therefore, the g-metric is invariant under 
the action of the (locally) timelike Killing vector <9/<9r, but the Stuckelberg field T and the 
f-metric are not invariant. As a result, the timelike isometry is not shared by both metrics. 

As the solution T = t is not static, it is unlikely to describe the “ground state” of 
the theory. This is probably the reason why this solution was found to be unstable [10]. 
Therefore, we need to consider other solutions for T. 

V. OTHER SOLUTIONS 

Solutions of the T-equation ( d t T) 2 — ( d r T ) 2 = 1 can be constructed in different ways. A 
fairly general solution containing an arbitrary function bL(£) is given by [16], 

T = cosh(£) t + sinh(£) r + W (£), 

dW(£) 

0 = sinh(£) t + cosh(£) r H-——, (5.1) 

where the second line implicitly determines the dependence of £ on t, r. Together with (4.2), 
this gives if not all but probably almost all solutions. However, this formula is difficult to 
use since one cannot explicitly determine £(f, r) for a generic W(£). 

The T-equation can also be integrated by applying the method of characteristics [16], 
which has a simple geometric interpretation. Let us consider the 2D Minkowski space 
spanned by x a = {t, r} with the metric g a b = diag[l, —1]. The T-equation reads g ab d a TdbT = 
(VT, VT) = 1. Let 7 = x a (s ) be a spacelike curve and T is constant along it. At every 
point of 7 there is a unit timelike normal n such that ( n,n) = 1 and (n,d/ds) = 0. The 
T-equation is equivalent to dT/dn = 1 [16]. 


IL + 

2.j 


(4.12) 
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FIG. 1. The T-equation in the method of characteristics. 


This allows one to pass from 7 where T = T( 7) to a neighboring curve 7 where T = T( 7) 
(see Fig.l) and so on, thereby extending T to the whole of the space. The solution is therefore 
defined, up to an additive constant, by the choice of the initial curve 7. For example, the 
solution (4.2) can be obtained by choosing 7 to be a straight line. 

In practice solutions of (dtT) 2 — (d r T) 2 = 1 can be obtained by changing the variables and 
then separating them [17]. Let us illustrate the method by passing to the static coordinates 
(4.11), in which case the T-equation becomes 


1 fdT\ 2 E fdT\ 2 

E V ~dr ) ~ 1 - E V ) ' 


(5.2) 


It is easy to see that T(r, p) given by (4.13) fulfills this equation, but now we can obtain also 
other solutions, in particular those for which dT does not depend on time and the f-metric is 
static. The most general solution of (5.2) of this type is obtained by separating the variables, 


T = \/l + g 2 7-+ J \/g 2 + p 2 , 


(5.3) 


where q is an integration constant. If q = 0 then the solution becomes especially simple, 

T = r + J ^ - p = V - p, (5.4) 

and choosing V and p as coordinates, the two metrics become 


ds 2 y = a 2 {-E dV 2 + 2dVdp + p 2 dQ 2 }, 
ds 2 = u 2 a 2 {—dV 2 + 2dVdp + p 2 dfl 2 }. 


(5.5) 
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VI. ENERGY 

One can compute the canonical energy for systems with non-trivial Stuckelberg fields in 
the same way as this is done in the unitary gauge [18]. The computation will be presented 
separately [19] but its result is as follows. For a solution expressed in the static coordinates 
t, p the energy on a hypersurface of constant r is 



with the radial energy density 

£ = u 2 P 2 (u)p 2 d r T. (6.7) 

Applying this formula to the T = t solution (4.13) gives the time-dependent value, 

£ = u 2 P 2 (u)p 2 \/l — p 2 cosh(r), (6.8) 

which indicates once again that this solution cannot describe the ground state. On the other 
hand, the energy will be time independent if d r T is time independent, but all such solutions 
are given by (5.3), in which case 

£ = u 2 P 2 {u)\/l + q 2 p 2 . (6.9) 

This corresponds to the constant volume energy density 

e = u 2 P 2 {u)y/l + q 2 , (6.10) 

and the total energy is E = eV where V is the (infinite) volume of the 3-space. We remember 
that u is a solution of the algebraic equation P\(u) = 0, therefore, depending on choice of u 
and also on values of the parameters &*,, the energy can be positive, negative, or zero. 

The actual value of the background energy is probably not so important, but it is impor¬ 
tant to know if the energy is minimal or not. We conjecture that the static solutions (5.3) 
correspond to the energy minima and are therefore stable. Therefore, they are candidates 
for describing the de Sitter ground state in the theory. To prove the conjecture will require 
to resolve the constraints and to compute the energy for deformations of the background 
[18]. We presently have partial results supporting our conjecture, but the detailed analysis 
will be presented elsewhere [19]. 
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VII. CONCLUSIONS 


We have shown that the de Sitter vacua in the dRGT theory are labeled by solutions of 
( d t T ) 2 — ( d r T ) 2 = 1. The simplest solution T = t is manifestly homogeneous and isotropic 


when written in the open chart, but it is unstable. Therefore, one should study other 
solutions. One could worry that other solutions will not be FLRW because their reference 
metric is inhomogeneous. However, as we have seen, this is not necessarily the case, as 
the reference metric can look inhomogeneous in some coordinates while sharing common 
translational isometries with the physical metric. 

The important issue is the number of common isometries of the two metrics. Since each 
of them describes a maximal symmetry space, each metric has ten isometries, some of which 
can be common, as for example the SO(3) rotational isometries. The number of common 
isometries depends on choice of T, for example for T — t this number is six, but the same 
can be true for other choices of T as well. 

Requiring the timelike isometry to be common for both metrics reduces the set of solutions 
to a one-parameter family (5.3). These solutions are physically distinguished since only for 
them the energy is time-independent. We conjecture that these solutions are stable and 
describe therefore the de Sitter ground state of the theory. The stability will follow if one 
shows that the energy increases for deformations of the de Sitter background, but such an 
analysis goes beyond the scope of the present paper and will be reported separately [19]. 
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Appendix A: T-equation in general coordinates 


The most general spherically symmetric g-metric can be represented as 



(A.l) 
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where N, (3, A, R are functions of two coordinates t, r. Choosing the f-metric to be 

dsj = a 1 2 3 4 5 6 7 8 u 2 {-dT 2 (t,r ) + dR 2 + R 2 dtt 2 } (A.2) 

and analyzing the field equation G^ = T^ u component by component, one finds [9] that 
they reduce to G^ + Ag, w = 0 with A = Pq{u) provided that Pi(u) = 0 and a 2 = 3/A, and 
if (p 2 + b 3 u) Y — 0 hence either b 2 + b 3 u = 0 or Y = 0. Here 

Y = (f — fiT' + NAR'^j 2 - (r - /3R’ + NAT'^j 2 - A(TR' - RT') + N^j 2 , (A.3) 

where the dot and the prime denote the derivatives with respect to t and r, respectively. 

If the g-metric is de Sitter and t , r coincide with the t, r coordinates of the ambient 5D 
Minkowski space used in (3.3), then one has 

N=-A=, R=VlTi r A - NR, = (A.4) 

Vl + 1 2 l + R 

Inserting this to the condition Y — 0 gives equation (3.20) for T(t,r) in the main text. The 
other possibility is to set b 2 + b 3 u = 0 [5],[20], which restricts the values of the parameters 
bk but the function T(t,r ) then remains arbitrary, which presumably indicates some hidden 
gauge invariance. 
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